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We investigate the growth index parameter 7 and the time variation of the gravitational constant 
G e s by using the currently available growth function f(z) data at different redshifts, with and 
without scaling to the fiducial ACDM model. We inquire the four different models of 7 including a 
constant 7. From a \ 2 minimization, we constrain the parameter spaces of models and show that 
ACDM model is excluded by 1-a level from current f(z) data. G e a is different from the Newton's 
gravitational constant Gn in modified gravity theories and interestingly, the current data shows 
that G ff 7^ Gn at 2 > 0.2 ~ 0.3 with 3-cr level. From these, we conclude that Einstein's General 
Relativity with ACDM is ruled out by 99 % confidence level from large scale structure observations. 

PACS numbers: 95.36.+X, 98.65.-r, 98.80.-k. 



1. INTRODUCTION 

Since the discovery of the current accelerating expansion of the Universe, both a dark energy (DE) and 
a modification of gravity (MG) at cosmological scales are the most commonly proposed candidates for 
explaining cosmic acceleration. It has been proposed that complementary observations of both geometrical 
tests and the cosmological evolution of the large scale structure (LSS) formation might reveal the origin of 
the cosmic acceleration among those models. 

It has been known that the flat ACDM cosmological model which is named to the so-called concordance 
model of Big Bang cosmology is consistent with the majority of current cosmological observations includ- 
ing both background evolution (e.g. supernovae type la, cosmic microwave background, baryon acoustic 
oscillation, etc) and LSS growth function (e.g. redshift space distortions, growth rate of clustering, weak 
lensing, etc). We re-examine the consistency of the ACDM model by using the current growth function 
data. We also investigate any time variation of the effective gravitational constant from the same data. 



2. CONSTRAINTS ON GROWTH FUNCTION 

The sub-horizon scales linear perturbation of the matter (<5 m = dp m / p m ) for many DE and MG models 
is governed by 

+ 2HS m = AttG 

cffPm^m ? 

(2.1) 

where dot means the derivatives with respect to cosmic time t, p m is the mean matter density, H is the 
Hubble parameter given by the Friedmann equation (H 2 = 8 g N p C r)i and G c s is the effective gravita- 
tional constant obtained from the Poisson equation. G e g becomes Newton's gravitational constant Gn in 
Einstein's general relativity. We change the variable in Eq. (|2.1j) from t to In a to obtain 
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where u is the effective equation of state of dark energy and VL m {a) — p m (a)/p c 
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The growth function, / is defined as the logarithmic derivative of the linear growth factor, D which is 
the growing mode solution of 5 m 
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where 7 is the so-called growth index parameter. For uj = —1/3 or —1 dark energy models, there exit the 
exact analytic solutions of 7 at present epoch [l| 
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where T is the gamma function, F is the hypergeometric function, and f2 m o is the present value of matter 

energy density contrast. Thus, for ACDM 7o = l n [-|r! m0 + 3fi| o r[^]/(r[|]F[§, §, f , 1 -0^])] /ln[fi m0 ] 

and varies from 0.558 to 0.553 for r2 m0 = (0.20, 0.35). There are also exact analytic solutions of D for general 
constant w dark energy models and thus one can obtain analytic form of 7 by using Eq. (£3D HI- 

In general, 7 can be a function of time and Eq. (|2.2[) can be rewritten by using / and 7 
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We adopt the currently available data of the growth function / = f3b derived from the redshift space 
distortion parameter fi(z) and the linear bias b(z) in order to constrain the free parameter 7. However, the 
data referred as "obs" in Table U have the different normalization (as) and present matter density (f2 m o). 
Thus, we adopt a rescaling method given in Q in order to transform different growth function to the same 
fiducial cosmology (ie., flat ACDM with f2 m o = 0.273 and as = 0.811) and denote the quantities obtained 
from this method with the subscription "res". If one assumes that both the redshift space distortion 
arameter /3 and the a$ are independent of fiducial model, then one can obtain the relation (see detail in 



/res(z) = /obs (2) 



O"8,obs -D b s (z) 

~V S D(z) 



(2.6) 



In Table HI we show the currently available values of growth function f bs(z) at the different redshifts 
which include 6 degree Field Galaxy Survey (6dFGS) at z — 0.067 [5(, 2 degree Field Galaxy Redshift 
Survey (2dFGRS) at z = 0.15, VIMOS-VLT Deep Survey (VVDS) at z = 0.77 % 2dFGRS at z = 0.17 
0, Sloan Digital Sky Survey (SDSS) luminous red galaxies (LRGs) at z = 0.35, 2dF-SDSS LRG and 
Quasi Stellar Objects (QSO) Survey (2SLAQ) at z = 0.55 9], 2dF QSO Redshift Survey (2QZ) and the 
2dF-SDSS LRG and QSO Survey (2SLAQ) at z = 1.40 ;10|, Large Sample of UVES QSO Absorption 
Spectra (LUQAS) at z = 2.42 [ll|, Lyman-a forest in the SDSS at z = 3.00 [II], and WiggleZ Dark Energy 
Survey at z = 0.22, 0.41, 0.44, 0.60, 0.73, 0.78 [T3, [3- We also derive the ratio of the effective gravitational 
constant to the Newton's one from / ODS and represent them as G e s/G-^\ bs- We also list the rescaled values 
of / res and Gcff/G N |r OS - 

We constraint parametrization of 7(2) with the data given in Table HI We use phenomenological 
parametrization of 7(2) by using the general relation 

l{ z ) = 7o + laP(z) . (2.7) 

We define Model 1 [ll|, 2 [lj|, 3, and 4 with p(z) = (z, — j-^, j-pj), and 0, respectively. Model 3 is also 
investigated, because 7 a is proportional to a in the constant ui models [l|. 
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TABLE I: The currently available data for the growth functions /obs and their re-scaled values / ros . Also we show 
the derived values of G e ff/GN|obs and G e ff /Glares from / b s and / rea , respectively. 



We perform a standard \ 2 minimization with two quantities / and g = Ggs/Gn to constrain the model 
parameters 70 and 7 a . 
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where N is the number of observations, O bs(thcor) is the observational (theoretical) values of observations, 
p is the fitted parameters, and Ui is the uncertainties of observations. We repeat y 2 minimization by using 
"observed" and "rescaled" values separately for different models. 

First, we constrain the 7 by using / h s with 13 (without asterisk) and 15 data points from Table |U We 
separate the two data points from others with asterisk, because the data obtained from the Lyman-a forest 
at z = 3.0 gives / > 1 which is physically unreasonable. Also, 2SLAQ data at z = 0.55 produces too large 
G e ff /Gn I obs- We define Data 1 and Data 2 for 13 and 15 data points, respectively. We define the reduced 
X 2 as Xrod = X 2 /dof with degree of freedom dof = (N — n — 1) where n is the number of fitted parameters. 
The results for \ 2 minimization of / Q b s are shown in Table ILTl and Fig. [TJ Ml provides (70,7a) = (0.47 ± 
0.07, 0.34 ± 0.20) with xLiJ d °f = 7.42/10 at 1-cr level for Data 1 and (70, 7a) = (0.48 ± 0.07, 0.32 ± 0.20) 
with Xmin/ d°f = 8.91/12 for Data 2, respectively. Data 1 gives the almost same good fit as Data 2. For M2, 
( 7Q) 7o ) = (0.45 ± 0.09, 0.57 ± 0.37) with x 2 m \Jdof = 7.85/10 and ( 7o , 7a ) = (0.46 ± 0.09, 0.55 ± 0.36) with 
xLiJdof = 9.22/12 for Data 1 and 2, respectively. If we use M3, then (70,7a) = (1.02±0.29, -0.57±0.37) 
with xLJdof = 7.85/10 and (70,7a) = (1.01 ± 0.29, -0.55 ± 0.36) with xLiJ d °f = 9-23/12 for Data 
1 and 2, respectively. M4 gives 70 = 0.58 ± 0.04 with xUJ do f = 10-21/10 and 7o = 0.58 ± 0.04 with 
Xmin/dof — 11.50/12 for Data 1 and 2, respectively. Thus, Ml is the best fit model and M4 is the worst 
one. M3 shows the almost same good fit as M2 but with large errors in 70. In Fig. [1] we show the 1 and 
2-cr likelihood contours of (70,7a) plane for Ml, M2, and M3 with Data 1 in the first row. We also show 
the same contours plots of same models with Data 2 in the second row. As shown in Fig. [TJ 7 values of 
all models exclude the ACDM expected value of 7(f2 m o = 0.273) = 0.555 at 1-cr level except M4. Also one 
of well known large extra dimension models of Dvali, Gabadadze, and Porrati (DGP) with 7 = 11/16 is 
excluded for all models at 2-cr level. However, 7 value of DGP is obtained from the different background 
evolution (i.e with oj different from -1) and thus one is not able to conclude for the validity of DGP model 
with data in Table HI 



Now, we repeat the constraint of 7 by using / res with 13 (without asterisk) and 15 data points from Table 
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TABLE II: 70 and j a obtained from \ 2 minimization of / b s for the different models with Data 1 and 2. 
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FIG. 1: a) Likelihood contours for l-a (solid lines) and 2-cr (dotted lines) confidence levels of the (70,7a) plane in 
the case of Model 1 (first column), Model 2 (second column), and Model 3 (third column). The first row is obtained 
with using Data 1. b) Same contours for the same models by using Data 2. 



HI We define Data 3 and Data 4 for 13 and 15 data points, respectively. The results for x 2 minimization of 
/ res are shown in Table |TTT] and Fig. [2 Ml provides (7o,7 Q ) = (0.50 ± 0.07, 0.34 ± 0.20) at l-a level with 
Xlajdof = 2.83/10 for Data 3 and ( 7o , 7„) = (0.50±0.07, 0.30±0.20) with Xmin/ rfo / = 5.80/12 for Data 4, 
respectively. As expected Data 3 gives the better fit than Data 4. For M2, (70, 7J = (0.48±0.08, 0.56±0.35) 
with x 2 m Jdof = 3.27/10 and (70,7a) = (0.49 ± 0.08, 0.51 ± 0.35) with x 2 m Jdof = 6.12/12 for Data 3 and 
4, respectively. If we use M3, then (70,7a) = (1-04 ± 0.28, -0.56 ± 0.35) with xLiJdof = 3.27/10 and 
( 70;7a ) = (1.00 ± 0.28, -0.51 ± 0.35) with xldjdof = 6.12/12 for Data 3 and 4, respectively. M4 gives 
7o = 0.60 ± 0.04 with x 2 min /dof = 5.73/10 and 7o = 0.59 ± 0.04 with X 2 m iJdof = 8.22/12 for Data 3 and 
4, respectively. Again Ml is the best fit model and M4 is the worst one. In Fig. [5J we show the 1 and 2-cr 
likelihood contours of (70,7a) plane for Ml, M2, and M3 with Data 3 in the first row. We also show the 
same contours plots of same models with Data 4 in the second row. As shown in Fig. [2l 7 values of all 
models exclude the ACDM expected value of 7(f2 m o — 0.273) = 0.555 at l-a level except M4. Also DGP 
with 7 = 11/16 is marginally excluded for all models at 2-a level. Compared to the constraints from Data 
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TABLE III: 70 and y a obtained from \ 2 minimization of / rea for the different models by using Data 3 and 4. 
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FIG. 2: a) Likelihood contours for 1-a (solid lines) and 2-a (dotted lines) confidence levels of the (70,7a) plane in 
the case of Model 1 (first column), Model 2 (second column), and Model 3 (third column). The first row is obtained 
with using Data 3. b) Same contours for the same models by using Data 4. 



1 and 2, / ros give the better fits for all models. The above results are quite similar to the recent analysis 
by using /cr 8 data fl7| . 



CONSTRAINTS ON TIME VARYING G 



Equation 



can also be used to check any time variation of the effective gravitational constant 
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TABLE IV: 70 and "f a obtained from \ 2 minimization of g ^ B for the different models with Data 1 and 2. 

where 7 can be obtained from the data by using its definition 7 = ln[/ b s (rcs)]/ m [^m.obs(rcs)] for the observed 
(rescaled) quantities . Thus, is obtained from Eq. (|3.1[) by adopting the fiducial ACDM model where 
we assume that 7 is a constant. The observed (rescaled) values of g bs(res) are listed in Table Q] We repeat 
the x 2 minimization to constrain the 7 parameters by using data from g. We separate Data 1 (3) and 2 
(4) with and without g b s (res) values at z — 0.55 and 3.00 again. 

We show the results for x 2 minimization of g b s in Table ITVl and Fig. [3J Ml provides (70, j a ) = (0.73 ± 
0.05, 0.82 ±0.22) at l-a level with xLiJdof = 27.65/10 for Data 1 and (70,7a) = (0.72 ± 0.05, 0.81 ± 0.22) 
with Xmin/dof = 28.67/12 for Data 2, respectively. Data 2 gives the better fit than Data 1. When we use 
M2, (70, 7a) = (0.70 ± 0.05, 0.67 ± 0.26) with xLJdof = 40.38/10 and (70, 7a) = (0.69 ± 0.05, 0.65 ± 0.26) 
with x'LiJdof = 41.25/12 for Data 1 Dnd 2, respectively. For M3, (70, la) = (0.96 ± 0.09, -0.65 ± 0.19) 
with x'LiJdof = 39.09/10 and (70, 7a) = (0.95 ± 0.09, -0.64 ± 0.19) with x'LiJdof = 39.88/12 for Data 1 
and 2, respectively. From M4, we obtain 70 = 0.72 ± 0.05 with xLiJdof = 48.13/10 and 70 = 0.72 ± 0.05 
with Xmin/dof = 48.70/12 for Data 1 and 2, respectively. Again, Ml is the best fit and M4 is the worst 
one. In Fig. [31 we show the l-a (inner darkest shaded regions), 2-er (dark shaded regions) and 3-er (light 
shaded ones) best fit form contours of g(z) for Ml, M2, and M3 with Data 1 in the first row. We also show 
the same contours plots of same models with Data 2 in the second row. Interestingly, G s/Gn deviates 
from ACDM expected value 1 from z > 0.2 ~ 0.3 for all models at 3-cr level. DGP expected value of g{z) 
is also shown as the dot-dashed line. However, g{z) of DGP is obtained from the effective equation of state 
of dark energy component oj = —1/(1 + O m (z)) which is different from the fiducial model and thus it is 
impossible to investigate the viability of DGP model with this data. 



Also, we repeat the x 2 minimization of <7 res and the results are shown in Table IVl and Fig. 0] Ml provides 
(70, 7 ) = (0.64 ±0.05, 0.61 ±0.23) with x'LiJdof = 4.46/10 for Data 3 and (70,7a) = (0.64 ± 0.05, 0.60 ± 
0.23) with xLiJdof = 5.41/12 for Data 4, respectively. For M2, we obtain (70,7a) = (0.62 ± 0.05, 0.68 ± 
0.26) with xLiJdof = 5.83/10 and (70,7a) = (0.62 ± 0.05, 0.67 ± 0.26) with xLiJdof = 6.77/12 for 
Data 3 and 4, respectively. For M3, (70,7a) = (0.91 ± 0.08, -0.67 ± 0.20) with xLiJdof = 5.81/10 and 
( 70;7a ) = (0.90 ± 0.08, -0.67 ± 0.20) with xLiJdof = 6.74/12 for Data 3 and 4, respectively. From M4, 
we obtain 7o = 0.69 ± 0.05 with xLiJdof = 13.20/10 and 7o = 0.69 ± 0.05 with xLiJdof = 14.05/12 
for Data 3 and 4, respectively. Again, Ml is the best fit and M4 is the worst one. In Fig. 01 we show 
the 1-er (inner darkest shaded regions), 2-er (dark shaded regions) and 3-er (light shaded ones) best fit form 
contours of g(z) for Ml, M2, and M3 with Data 3 in the first row. We also show the same contours plots of 
same models with Data 4 in the second row. Again, G e s / Gn deviates from ACDM expected value 1 from 
z > 0.3 for all models at 3-a. 



4. CONCLUSIONS 

After we have compiled a currently available dataset of the growth function from the fiducial concordance 
ACDM model for the background evolution, we also re-scale the values of the observed growth function. We 
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FIG. 3: a) The best fit form of Geff/Gr-jlobs for l-cr (darker shaded regions), 2-cr (dark shaded regions), and 3-cr 
(light shaded regions) errors in the case of Model 1 (first column), Model 2 (second column), and Model 3 (third 
column). The first row is obtained with using Data 1. b) Same best fit form of G e s/G-N the same models with using 
Data 2. 
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Xrcd.min 


Ml 


0.64 ± 0.05 


0.61 ± 0.23 


4.46 


0.45 


0.64 ± 0.05 


0.60 ± 0.23 


5.41 


0.45 


M2 


0.62 ± 0.05 


0.68 ± 0.26 


5.83 


0.58 


0.62 ± 0.05 


0.67 ± 0.26 


6.77 


0.56 


M3 


0.91 ± 0.08 


-0.67 ± 0.20 


5.81 


0.58 


0.90 ± 0.08 


-0.67 ± 0.20 


6.74 


0.56 


M4 


0.69 ± 0.05 





13.20 


1.32 


0.69 ± 0.05 





14.05 


1.17 



TABLE V: 70 and 7 a obtained from \ 2 minimization of g les for the different models with Data 3 and 4. 



constrain the several phenomenologically parameterized growth index parameter 7 allowing time variation 
of it. Theoretically, almost the constant 7 is expected from the ACDM model, but the best model from 
original data produces 7(2) = 70 + J a z with the best fit values (70,7a) = (0.48 ± 0.07,0.32 ± 0.20) and 
(0.47 ± 0.07, 0.34 ± 0.20) with and without 2 SLAQ and Lyman-a data, respectively. The reduced \ 2 
values are both 0.74. When we use the rescaled values of /, then the best fit values are (70,7a) = 
(0.50 ± 0.07,0.30 ± 0.20) and (0.50 ± 0.07,0.34 ± 0.20) with and without 2 SLAQ and Lyman-a data, 
respectively. The reduced x 2 values are 0.48 and 0.28, respectively. Even though the constant 7 model 
includes the ACDM expected value of 7A(f2 m o = 0.273) = 0.555 at l-cr level, this gives the worse fit 
compared to other models for all data. 

We also investigate the time variation of the effective gravitational constant proposed by modified gravity 
theories. It is interesting to find that the effective gravitational constant is smaller than the Newton's 
gravitational constant Gn for z > 0.2 ~ 0.3 for all considered models at 3-cr level. Thus, we can conclude 
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z z z 

FIG. 4: a) The best fit form of G e ff/GN for 1-a (darker shaded regions), 2-a (dark shaded regions), and 3-er (light 
shaded regions) errors in the case of Model 1 (first column), Model 2 (second column), and Model 3 (third column). 
The first row is obtained with using Data 3. b) Same best fit form of G e ff/GN the same models with using Data 4. 



that the current large scale structure formation shows the deviation from Einstein's General Relativity 
ACDM with 99 % confidence level. 

The current growth rate observations show the deviation from Einstein's general relativity. By using 
the previous data with 20 — 40 % accuracies and WiggleZ galaxy survey within 9 — 17 %, the errors in 7 
parameters are about 15 % and 60 % for 70 and j a , respectively. 
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